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Abstract 

The properties of principal null directions of a perturbed black hole are 
investigated. It shown that principal null directions are directly observ- 
able quantities characterizing the space-time. A definition of a perturbed 
space-time, generalizing that given by Stewart and Walker is proposed. 
This more general framework allows one to include descriptions of a given 
space-time other than by a pair (M, g) where M is a four-dimensional dif- 
ferential manifold and g a Lorentz metric. Examples of alternative char- 
acterizations are the curvature representation of Karlhede and others, the 
Newman-Penrose representation or observable quantities involving prin- 
cipal null directions. The conditions are studied under which the various 
alternative choices of observables provide equivalent descriptions of the 
space-time. 



1 Introduction 

The necessary and sufficient condition for a space-time to be algebraically special 
in the Petrov classification is P| that the curvature invariants I and J satisfy 

I 3 - 27J 2 = 0. (1) 

The curvature invariants are the determinantal expressions in the components 
of the Weyl spinor [2] 
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A Newman-Penrose (NP) notation for the components of a symmetric spinor is 
used here, i.e., 

*m = *0...1 (3) 



with m the number of spinor indices 1. 

Baker and Campanclli 3 introduce the speciality index of the space-time as 
the ratio 



For a Kerr-Newman black hole in the Kinnersley tetrad gy, the only nonvanish- 
ing component of the Weyl spinor is ^2- Expressions J2J) then take simple forms 
and the condition is satisfied (as it should for a type D space-time). Baker 
and Campanelli go on and expand the speciality index in powers of an arbitrary 
perturbation parameter e of the black hole: 



This is a perplexing result since it implies that perturbed black holes, in the 
first-order approximation, are algebraically special. There exists an extensive 
literature El 13 E] of 'algebraically special perturbations'. What are these 
then? It is our purpose in the present paper to shed some light on this apparent 
controversy by examining the properties of the principal null directions. The 
clarification of this point leads one to fundamental issues such as the notion of 
the perturbation of a space-time. 

A general framework for perturbations of space-times has been discussed by 
Stewart and Walker [§]. They propose the following 

Definition 1 (Stewart- Walker) A space-time (M.',g') consisting of the 
manifold M 1 and metric g' is a perturbation of some given space-time (M,g) 
if there exists a smooth one-parameter family of space-times (M,g(X)) in a 
Hausdorff five- dimensional manifold M^ 5 -* which contains (A4,g) = (A4,g(0)) 
and (M' , g') = (A4, g (e)) . The unperturbed or background space-time (M,g) 
is given by A = 0. In the parameter interval A G (— e,e), the g ab (A) define 

a piecewise smooth tensor field g a @ on of signature (OH ). The 

singular hypersurfaces g a @ (d\)g = are then the original space-time manifolds 



A weakness of this definition is that it is too restrictive to be applicable to 
some important works in perturbation theory. Many treatments of black-hole 
perturbations characterize the space-time by using curvature components and 
other quantities which are not considered in the definition. As an example, 
Chandrasekhar JO] characterizes a perturbed black hole such that the dyad 
components of the Weyl spinor plus the optical scalars differ only by first-order 
terms from their form in the background space-time. Such treatments are admis- 
sible only in a general framework in which the characterization of the space-time 
in more than one ways is allowed. For example, in the curvature representation 
jll) . a space-time is locally completely determined by the Riemann tensor and 
a finite number of its derivatives in a moving frame. 

It thus appears necessary to formulate the perturbation problem in a more 
general way than by the Stewart- Walker definition. If one wants to encompass 
the various alternative characterizations of the space-time, one may replace the 
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pair (M',g') in the definition by a pair (A4' , O') where O' is some complete set 
of measurable quantities characterizing the space-time. The complete charac- 
terization of the geometry is the central issue in what is known the equivalence 
problem Q31H1- 

Definition 2 A set of observables O on some open set U is said to locally 
completely characterize the space-time containing U if the differ entiable and 
metric structures are fixed uniquely on hi by the values of O. 

Obvious examples of the choice of O' are the perturbed metric g' or the 
quantities used in the curvature representation [TT]. An incomplete set of ob- 
servables, for example, in a space-time containing a number of different matter 
fields is obtained when dropping the subset of observables such as the field 
stresses of one of the fields from among the observables. 

The proposed new definition of a perturbed space-time raises some ques- 
tions which did not complicate the old definition. An important such question 
is whether two different choices of O' yield equivalent treatments of the pertur- 
bation problem. In the context of linear perturbations, we may formulate this 
as the requirement that the two choices O and O' of the observables are related 
by O' = M O where M is a nonsingular matrix. 

Our goal in this paper is to investigate the sets of observable quantities 
characterizing the space-time in order to clarify if different choices of these 
sets of observables will yield equivalent descriptions of the perturbations. In 
particular, we want to clarify how the various choices of observable quantities 
behave under perturbations of the system. As a first step on this route, in the 
next section we recapitulate the behavior of rotational metric perturbations in 
Hartle's theory |13] . In Sec. we briefly recall the concept of principal null 
directions of a symmetric spinor. The application of this general theory to the 
electromagnetic and gravitational fields in a Kerr-Newman black hole is reviewed 
in Sec. The behavior of principal null directions of the electromagnetic vs. 
gravitational field to first order in the perturbation parameter is worked out in 
Sees. Eland El respectively. 

We find, as an unexpected result of this investigation, that there are certain 
choices of the sets of observables which provide inequivalent descriptions of the 
space-time. The existence of these inequivalent sets is shown to have impor- 
tant consequences for the theory of gravitational radiation. This finding, when 
combined with the properties of principal null directions in the perturbative 
picture, makes it possible to reach the correct interpretation of 'algebraically 
special perturbations'. 

2 Illustrative example: rotational perturbations 

In the Hartle theory |13| of rotational perturbations, the metric of the un- 
perturbed vacuum is given by the Schwarzschild solution with mass m. The 



3 



perturbed space-time has the metric 

o , ( 2m\ 9 „ / 2m\ 1 „ 
ds 2 = (1 + 2h) ( 1 J di 2 - (1 + 2j) f 1 J dr 2 

- (1 + 2k) r 2 \dtf 2 + sin 2 & (dip - ujdtf] . (6) 

Here the perturbation functions h, j, k and u> are assumed to depend smoothly 
on the perturbation parameter, i.e., the angular velocity Q. Additional infor- 
mation on the form of these quantities is obtained from the discrete symmetries 
of the rotating system. Taking into account that the simultaneous reversal of 
time direction and the reversal of the sense of rotation is an exact symmetry, 
we have that the series expansion of the function to contains only odd powers 
of the angular velocity. The expansion of all other perturbation functions has 
only even powers in Q. Hence one finds that the only unknown function in the 
perturbation problem to first order in £1 is u>. In fact, the field equations yield 
an uncoupled, second-order linear differential equation for the function lu. Given 
the solution of this equation, one may proceed to solve the perturbation prob- 
lem to second order in the angular velocity f2. In the second approximation, 
the field equations are coupled second-order linear differential equations for the 
remaining perturbation functions. The function lu contributes quadratic source 
terms in these equations. 

The rotational perturbations of the Schwarzschild black hole illustrate our 
point: a description which takes into account the full physical properties of the 
perturbed system can yield restrictions on the form of the series expansion in 
the perturbation parameter. 



3 Principal null directions 

The principal spinors of a k- index symmetric spinor lpabc.k are defined as the 
nonzero spinors o A satisfying the condition 

Vabc...ko a o b o c ...o k = 0. (7) 

Let o° be the nonzero component and let us introduce the complex ratio 



o" 

The roots of the complex algebraic equation 



ip + kcp lZ + ( 2 ) ^ + ( 3 J ^ + •- + Vkzk = ( 9 ) 

[where we are using the notation for the components of the spinor] define the 
flagpoles o A o A . The principal null directions of <Pabc...k are represented, 
up to real multiplying factors, by these flagpoles. By the fundamental theorem 
of the algebra, there exist k roots. 



(8) 
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Figure 1: The angle subtended by two principal null directions for an observer 
is given by the line element ds of his celestial sphere. 



The various coincidences of the principal null directions are classified by the 
partitions of the number k. In the generic type, there are k distinct principal 
directions, represented by k vectors of the light cone in the tangent space T p 
at point p of the space-time (Figure 1). All other types are called algebraically 
special. 

Principal null directions at any given point are directly observable quantities 
of the space-time. An observer with world line through the point p can detect 
principal null directions by looking at a distant sphere [15] . In the lowest- 
order approximation, the observed image of the sphere is distorted to an ellipse. 
Kristian and Sachs ^H] characterize the distortion by the ratio of the major and 
minor axes of the ellipse: 



depending both on the distance R and the direction <p) of the source. To the 
present accuracy, the distance R can be equivalently taken to be the luminosity 
distance, the distance by the apparent size or area distance. The coefficient 
C (R, <p) is the projection of the Weyl tensor in the tangent plane of the 
celestial sphere |lfi) . At the principal null directions, the coefficient C of the 
distortion has the limiting value liniR = o C (i?, i?, tp) = 0. 

The angle subtended by a pair of these directions is given as the invariant 
length s on the unit sphere S 2 with metric 



D = 



hc{R,$,<p) R 2 



(10) 



ds 2 = dti 2 +sin 2 i?<V. 



(11) 



In an orthogonal frame {t, x, y, z}, the position of point p of Minkowski space- 
time is Tt + Xx+Yy + Zz. The sphere S 2 is represented ^3] by the intersection 



of the past null cone with the hyperplane T = — 1. The polar coordinates i? 
and ip are related to the complex stereographic coordinates z and ~z by the 
transformation 

z = e iV cot^. (12) 
The metric of S 2 in stereographic coordinates has the form 

ds 2 = ^rdzdz. (13) 

(l + zzf 

The coordinate system (z,z) is singular at the north pole where z — > oo, and 
the coordinate distortion due to the presence of the metric coefficient P = 
2T72 (1 + zz) increases indefinitely as one approaches the pole. 

There exists a sizable literature of the properties of principal null directions. 
Penrose and Rindler devote a chapter of their monograph to a detailed study 
of this subject. Gunnarsen et al. |15| develop a numerical approach, based on 
the d'Inverno-Russel Clark version of the Ferrari formula, for computing the 
gravitational principal null directions. They employ this method for the Kastor- 
Traschen metrics containing colliding black holes. The status of observations 
of principal null directions, along with various proposed techniques has been 
reviewed by Chrobok and Perlick 1 7 . 



4 The Kerr-Newman black hole 

The Kerr-Newman black hole with mass m, rotation parameter a and electric 
charge e has the metric 

ds 2 =(l- 2m ^_ e ^ (dt -asm 2 tfd(p) 2 

+ 2 (dt - a sin 2 <&d<p) (dr + a sin 2 <&d(p) (14) 
- CC (d$ 2 + sin 2 tfdtp 2 ) 

with the four-potential 

A = — =(dt — a sin 2 'ddip^j (15) 

where 

C = r — iacosi?. (16) 

The coordinates t and r run from —00 to 00 while $ and if are coordinates on 
a 2-sphere such that tp is periodic with period 2n and $ ranges from to n. In 
the generic case, the expression 

A = r 2 + a 2 - 2mr + e 2 (17) 

has two distinct roots. 
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The two double gravitational principal null directions in the Kinnersley 
tetrad are at the south versus the north poles of S 2 , i.e., at the values of 
the stereographic coordinates z = and z = oo. Since the coordinate system is 
singular at the north pole, it is more advantageous to use the null tetrad [2] 



D = {*— = — 

dx a dr 
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In the NP notation, the Maxwell tensor components are 
$ = F ab e a m b = 



$! = ^F ab (£ a rr 
$ 2 = F ab m a n b = 



2 1 / 2 C 2 



5 sin $ 



C 3 



The significant components of the Weyl curvature C a bcd are 

" 2 — m£ 



*2 
*3 



C 3 C 



-3ia sin?? 



m£ 



2 



2 • 2.a m C- e 



* 4 = 3a 2 sin 2 § 

and the remaining two components vanish, 
*o = -C abcd e a m b £ c m d = , tt a 



C 5 C 



-c abcd e a n b e c m d 







(18) 



(19) 



(20) 



The curvature singularities are at r = and i? = 7r/2 . 

Let the complex ratio of the Maxwell tensor components be denoted 



J>2_ 

2$i 



zasint? 
2 1 / 2 Q ' 



In this notation, the Weyl tensor components satisfy 



(21) 



(22) 



The double gravitational principal null directions are at the finite coordinate 
values z = and z = —1/w. 
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The observation of principal null directions in a Kerr-Newman black hole can 
be accomplished in an alternative way |18j to the one described in the previous 
section. There exist special null congruences ('photon trajectories') consisting 
of the integral curves of principal null directions, characterized by vanishing first 
curvature k and shear a. For these curves, the polar angle $ is constant. The 
axis value of the angular momentum is L z = aE sin 2 i?, with E the conserved 
energy of the motion, and the Carter separation constant vanishes, K = 0. 



5 Electromagnetic perturbations 

For later reference, we consider here the perturbations of the electromagnetic 
field in the neighborhood of a black hole. Our first choice for the (incomplete) 
set of observable quantities is the Maxwell field components. The unperturbed 
field is given by the components (|19|l . The perturbed Maxwell field has the 
expansion 

$ = e$£ ] + O (e 2 ) 

*i = $i 0) + + O (e 2 ) (23) 
4> 2 = <4 0) + e*^ + O (s 2 ) 

where the parenthesized superscripts indicate the degree in the perturbation 
parameter e. 

The principal null directions of the perturbed Maxwell field are given by the 
solutions of the quadratic equation 

z 2 + a 1 z + a = (24) 

where the coefficients are defined 

A possible alternative choice of the observable quantities is the set of principal 
null directions. We now want to examine if these two choices of the sets of 
observables are equivalent. The relation between the coefficients ai and the 
roots z\ and Z2 is given by the Viete formulae 

= z 1 + z 2 (26) 

Expanding all observable quantities in power series of e , we have 

a = ea^ + O (e 2 ) 

oi = af ) + ea/p + O (e 2 ) 

Zl =ez[ 1) +0(e 2 ) (27) 
z 2 = 4 (S >+e4 1) + 0(e t ) 
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where the unperturbed principal null directions are at z^ = and z^ = —1/w. 
Inserting in Eqs. Ij26(l . the unperturbed terms cancel and the first-order parts 
yield 

aW=_(^)+,W) (28) 

From these relations we see that the set (z^\ Z^j of first-order observables 

is linearly equivalent to the set ^a^\a^^j . Our conclusion is that the field 
stresses of the perturbed Maxwell tensor yield a description equivalent with the 
perturbed principal null directions. 



6 Gravitational perturbations 

In the literature of perturbed black holes ^U], the sets of observable quantities 
characterizing the gravitational field are selected from several available options. 
The most frequent two choices are the dyad components of the Weyl spinor in 
the NP formalism, and the components of the perturbed metrics. These two 
descriptors of the state are related to each other by the second-order differential 
equations embodying the field equations and the gauge conditions. Rather than 
investigating these complicated relations, we want to compare here the two 
choices of the observables given by the dyad components of the Weyl spinor and 
by the principal null directions. These alternatives are connected by the quartic 
equation 

*42 4 + 4^32 3 + 6*2^ 2 + 4*iz + *o = 0. (29) 

For the unperturbed black hole, the Weyl spinor components have the form 
120(1 and the four solutions for the principal null directions pairwise coincide: 
z^ = (0,0, —1/w, —1/w). We introduce the normalized coefficients 

for i = 0, ...,4. As in the previous section, we expand all observables of the 
perturbed space-time in powers of the perturbation parameter e, 

a k = 4 0) + ea ( k 1] + O (e 2 ) far k = 0, 3. (31) 

The two sets of observable quantities are now related by the symmetric expres- 
sions 

4 4 

-a 3 ^^Zi, a 2 = ^ ZiZk 

i—l i,k—l. i<k 

4 4 

i,k,l=l, i<k<l i—l 
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Inserting here the series expansions of the roots 

z 1 =ezf ) +0(e 2 ) 
z 2 = ez { 2 1] + O (e 2 ) 

z 3 = --+ezP +0(e 2 ) (33) 

z i = ~-+ezi 1] +0(e 2 ) 
w 

we get the linear relations for the first-order observables 





= 




(34) 


«« 


- — (z {1) 




(35) 








(36) 


4 1} 


= - 4 




(37) 



Hence we see that the choice of the principal null directions as the set of ob- 
servables is inequivalent to choosing the set of Weyl spinor components. By Eq. 
Ij34(l . for space-times with perturbed principal null directions, the coefficient ao 
is perturbed only at the second order in the parameter s. The second-order 
contribution, however, is nonvanishing and has the form = ■^tZ^Zq. 

Here we have an example of a situation when two choices of the sets of observ- 
able quantities yield inequivalent characterizations of the perturbed space-time. 
We then need to decide which of these should be used in the definition of the 
perturbed space-time. A significant portion of the existing literature uses the 
Weyl tensor components for characterization, with no regard to the behavior of 
the alternative set of observables such as the principal null directions. We argue, 
however, that two space-times having finite (or of lower-order) differences in 
some of their observable quantities should not be considered as exhibiting small 
perturbations. In the above situation, the two black hole space-times that have 
first-order differences of the curvature component "So do not satisfy the linear 
relation (|34|l . hence they are inequivalent in the sense that they must have their 
corresponding principal null directions pointing at finite angles. 

The correct interpretation of the above situation can be inferred from the 
behavior of rotational perturbations. When one choice of the set of observables 
gives rise to space-times which are characterized by finite differences in certain 
observables, then this indicates that pairs of corresponding quantities with in- 
finitesimal differences will only change in a higher-order approximation with the 
appropriate choice of the perturbation parameter. For rotational perturbations, 
one might carelessly expand the diagonal components of the metric with respect 
to some parameter (such as the square of the angular velocity) such that the 
perturbed state has first-order contributions in the diagonal metric. That this 
is not the right choice of the parameter, becomes clear when considering the 
behavior under reflections. 

10 



7 Gauge choices 



The foregoing discussion did not extend to the problem of choosing the tetrad 
gauge. In the exact description of a space-time, the choice of a null tetrad 
vector along a principal null direction leaves only a small discrete group of 
gauge symmetries permuting the principal directions. In the present section we 
show that the gauge group in the linear approximation is much larger, and in 
fact it has an important geometrical content. 

The perturbed tetrad is chosen such that the spinor o A is one of the four 
principal spinors of the Weyl curvature: 

*o = 0. (38) 

For the unperturbed black hole, two independent choices for the spinor o A are 
possible, corresponding to the two distinct double principal directions. 

In perturbation theory, the solutions of l|38[) remain, to some extent, unde- 
termined in the linear approximation. Rather than having a finite number of 
exact solutions of Eq. I|38[l i there exist infinitely many solutions of the problem. 
This can be seen from the expansions i|35|) of the symmetric expressions. Given 
the coefficients c£ , these equations have an infinite number of solutions for the 
roots 2- . The multitude of the principal null directions in the linear approx- 
imations gives rise to a group of dyad transformations which leaves condition 
unchanged. 

The quantity transforms under the infinitesimal dyad transformation 

o A ^o A + bi A , l a l a , (39) 

as follows: 

*o *o + 46*i. (40) 

Here b is an arbitrary but small complex multiplier function such that higher 
powers of b are negligible. Since the curvature quantity itself is small, the 
spinor o A remains a principal spinor of the curvature under the transformations 

Gauge symmetries arc normally considered as mathematical properties of 
our description. A gauge transformation is expected to alter such quantities 
of the description as the tetrad frame or a field potential while leaving the 
physically measurable quantities (the curvature or the field strength) invariant. 
From this point of view, the above setup where we choose a tetrad along a 
principal null direction is somewhat unusual. Here a gauge transformation of 
the type (|39|l amounts to selecting a congruence of principal null curves. The 
particular congruence thus selected is that with a tangent spinor o A . Thus a 
'gauge change' of the form l|39[l amounts to picking a new family of principal 
null curves. These families of curves may have different propagation properties. 
Hence the choice of the 'gauge' has significant consequences for the resulting 
behavior in the linear approximation. 

One of the earliest contributions to the literature of algebraically special per- 
turbations by Couch and Newman [S] takes advantage of an especially simple 
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form of the metric, and finds the perturbations of a Schwarzschild black hole. 
The linear algebraically special perturbations of the Kerr black hole are consid- 
ered by Wald |S] in the NP approach. Wald uses a gauge in which = ^1 = 0. 
For this particular pick of the principal null congruence, both the shear and the 
first curvature vanish. Wald's aim is to prove that also ^4 vanishes, in order 
to establish that the knowledge of ^0 alone suffices for a full description of the 
perturbed space-time. He fails to achieve this goal for modes with certain 'alge- 
braically special' frequencies. Note that one should proceed with caution when 
using the results of this paper. For example, on p. 1457, one finds this claim: 
'i.e., we obtain the linearized version of the Goldberg-Sachs theorem'. Now the 
Goldberg-Sachs theorem has been shown not to hold in the linearized theory 

m 

Chandrasekhar further investigates the properties of algebraically special 
perturbations [1] ■ He concludes his work with several, apparently technical, but 
at any rate, unanswered questions. These concern the relation of the Starobin- 
sky operator to complex conjugation, and to the special forms of the potential 
barriers occurring. These questions (he writes) are sheated in enigmas. A de- 
tailed investigation of the algebraically special frequencies is to be found in . 

Despite the deceptive simplicity achieved by the choice of the principal null 
congruence with = 0, this approach neither offers a complete description of 
the space-time nor can it be easily generalized for the presence of charge. 

8 Second-order perturbations 

In Section 6, we have shown that the component ^0 of the curvature does not 
receive a contribution in the linear approximation. In the light of existing earlier 
literature of black hole perurbations, it is an important question what are the 
equations governing the first nonvanishing contribution to the field ^q. 

Our goal in this section is to obtain an uncoupled differential equation for 
the second-order function \l/o- To this end, we express the derivatives 5p, 5k, 
Dt, Df3, Dp and Da from the NP equations (4.2k), (4.2b), (4.2c), (4.2e), (4.2h) 
and from the complex conjugate of (4. 2d), respectively. Similarly, we express 
D^2, 5^2; D^i and 5^>\ from the NP electrovacuum Bianchi identities (A3). 
Where necessary, second derivatives are untangled by use of the commutator 
[D, 6}. Acting with the commutator [D,S] on^!, the resulting equation has the 
form 

[55 -DA- Aa5 + 2aS + (4 7 - p) D + (p + 4p) A 

-8aa + 4£>7 - 46a - Apj - 16-yp + 4pp + 2* 2 ] *o (41) 

= -4(5- 6a) (cr#i) + 4 (A - 7 + p - 57) - 10*? . 

To second order, the terms in the square bracket on the left can be taken to 
have their value in the Kerr-Newman metric since they act on the second-order 
quantity *o . The terms on the right contain either of the second-order quantities 
ct^i, k^i or The operators acting on these second-order quantities assume 
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their values at the Kerr-Newman metric: 

V2 
' C 
m(C + C) 



□2*0 



iasin?? 



0_ 
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cos 2 $ 



sin 1? 



2a 



1 2 sin 2 tf^r - 2 



1 



sin 2 
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2is cosfl 



d_ 
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(r 2 + a 2 )— + (s + 2)r + 2 a cost? 
dr 



d_ 
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(42) 



K - > * 



(43) 



is the wave operator introduced in |2()| . 

Given the first-order solution of the perturbation problem, we have an uncou- 
pled linear differential equation for the unknown function ^o- The right-hand 
side is fully known and is to be treated as a source term. 

In the second approximation, the issue of the gauge choice arises again. This 
has been investigated in |23 and |22| . 



9 Conclusions 

We have argued above that perturbations of a black hole can be described in 
a sense which conforms to generally accepted criteria if the lowest-order contri- 
butions to the curvature component \I/o are of second order in the perturbation 
parameter. First-order contributions to Wo would result in finite changes in the 
principal null directions, which are directly accessible to observations. Before 
these disturbing results can be considered as fully consolidated, there is a clear 
need to clarify a number of related issues. We conclude this work with raising 
just one such example. 

The curvature component W4 remains invariant under transformations (|39|) . 
Under infinitesimal dyad transformations of the form 

. A . A , A „A . A 



l a + co a , o A ->o A , (44) 



with c an arbitrary small complex function, W4 transforms 

* 4 -> *4 + 4c* 3 . (45) 

In the Kinnersley tetrad, the spinor l a is chosen, as is o A , to be a principal spinor 
for the Kerr metric and both ^3 and W4 vanish. Repeating the argument 
with the symmetrical expressions (|35|l for W4 and projective coordinate z' = 
— 1/z, at first sight it would seem to be possible to choose the null tetrad for 
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the perturbed space-time such that the only non- vanishing tetrad component of 
the Weyl spinor is again \&2- This would imply that the perturbed space-time 
is again the Kerr black hole with trivial parameter changes. 

To find ways out from controversies like the one just described, either one 
has to give up insisting on the observability of principal null directions, or else 
to carry out a careful investigation of perturbations in a neighborhood of the 
coordinate singularity on the two-sphere. We do not see any ground, however, 
on which the observability of principal null directions should be given up. 
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